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Among the widely studied superconducting iron-pnictide compounds belonging to the Lnllll family 
(with Ln a lanthanide), a systematic investigation of the crossover region between the superconducting 
and the antiferromagnetic phase for the Ln = Nd case has been missing. We fill this gap by focusing on the 
intermediate doping regime of NdFeAsOj^F* by means of dc-magnetometry and muon-spin spectroscopy 
fiuSR) measurements. The long-range order we detect at low fluorine doping is replaced by short-range 
magnetic interactions at x = 0.08, where also superconductivity appears. In this case, longitudinal-field 
juSR experiments show clear evidence of slow magnetic fluctuations that disappear at low temperatures. 

This fluctuating component is ascribed to the glassy-like character of the magnetically ordered phase of 
NdFeAsO at intermediate fluorine doping. 

PACS numbers: 74.25.Dw, 74.25.Ha, 76.75.+i 


I. INTRODUCTION 

The simultaneous presence of superconductivity (SC) 
and magnetism (M), competing and/or coexisting in var¬ 
ious families of iron-based superconductors, has captured 
the interest of the condensed matter community for many 
years now. Particular attention has been devoted to the 
Lnllll family (Ln being a lanthanide), where the static 
long-range antiferromagnetic order (AFM), occurring in 
the FeAs planes, is progressively suppressed in favor of 
the superconducting state by carrier-doping of the parent 
compound. In the Lnllll family the M-SC crossover re¬ 
gion is characterized by the coexistence of a short-range 
AFM order and superconductivity on a nanoscopic length 
scale. This distinctive feature has been shown to be in¬ 
dependent of the doping agent, since it occurs both in FI- 
doped Lallll,- as well as on F-doped Ce-1111,— or Sm- 
Illli2 While this suggests that H or F act merely as elec¬ 
tron donors (although with different efficiencies), it leaves 
open the question of the lanthanide ions, i.e.: are the M-SC 
ground-state properties of the Lnl 111 family also indepen¬ 
dent of the Ln ion? 

To address this issue and to better understand the role 
of LnO planes on the occurrence of a M/SC ground state, 
we considered the case of NdFeAsO,_ r F A ., a family about 
which relatively fewer studies have been published inso¬ 
far, none of which addressing this important point. In 
particular, since the discovery of superconductivity in a 
Ndllll compound in 2008,- most of the studies focused 
on the evolution of the structural and superconducting 
properties induced either by F-doping,— or by simple oxy¬ 
gen deficiency^ The magnetic properties of the undoped 
compound and those of an x(F) = 0.06 under-doped sam¬ 
ple were reported in Refs. @ and [§, respectively. In the 
latter case it was emphasized that the dc magnetization re¬ 


sponse is dominated by the Curie-Weiss paramagnetism of 
the Nd 3+ ions, whereas no signs of long-range antiferro¬ 
magnetic order, due either to Nd 3+ or to Fe 2+ ions, could 
be detected below T c . Later on, Mossbauer-effect studies 
confirmed the absence of magnetically-ordered phases in 
an x(F) = 0.12 superconducting sampled /iSR measure¬ 
ments on either undoped,— or on F-dopedi2 samples were 
performed, too. In the former case they could detect an 
AFM long-range magnetic order with T N = 135 K and a 
local field B^ = 0.17T, while in the latter a gapped-like 
superfluid density was inferred. Unfortunately, due to the 
lack of good-quality, finely-tuned F-doped samples, a sys¬ 
tematic study of the most interesting region, that of the 
M-SC crossover, has been missing to date. 

In the present paper we address exactly this issue by in¬ 
vestigating via /iSR four different F-doped polycrystalline 
samples with x(F) = 0.03, 0.05, 0.08, and 0.25. In the first 
two cases we find a long-range AF order, that gradually ex¬ 
tinguishes as the F content increases. On the other hand, 
in the x(F) = 0.25 sample, we do not find any magneti¬ 
cally ordered phases. The intermediate x(F) = 0.08 case, 
instead, is particularly remarkable since here we detect 
the simultaneous presence of bulk superconductivity (be¬ 
low T c = 20 K) and of short-range magnetic order (below 
T n = 40 K). Interestingly, in an intermediate T c < T < T N 
temperature range, sizable spin-fluctuations with a slow 
dynamics dominate. We ascribe these magnetic spin fluc¬ 
tuations to the glassy character of the ordered magnetic 
phase. 

In a broader context, the current study of the Ndllll fam¬ 
ily fills a gap about the role played by the different rare 
earths in determining the properties of Lnllll compounds, 
therefore allowing detailed comparisons to be made. 
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II. EXPERIMENTAL 

A. Sample synthesis and structural characterization 

A series of NdFeAsO,_ x F x polycrystalline samples, with 
nominal fluorine content ranging from x = 0.03 to x = 
0.25, was prepared at ETH Zurich using cubic-anvil, high- 
pressure, and high-temperature techniques^ - — Due to the 
toxicity of arsenic, all procedures related to the sample 
preparation were performed in a glove box. Pellets contain¬ 
ing the high-purity (> 99.95%) precursors (NdAs, FeAs, 
Fe 2 0 3 , Fe, and NdF 3 ) were enclosed in a boron nitride con¬ 
tainer and placed inside a graphite heater. Six tungsten 
carbide anvils generated the required high pressure (3 GPa) 
on the assembly. While the pressure was kept constant, the 
temperature was ramped up in 2 h, from 20° C to the max¬ 
imum value of ~ 1350°C, and maintained there for 15 h. 
Following the temperature quenching, the pressure was re¬ 
leased and the sample removed. 

The structural characterization was performed by means 
of standard powder x-ray diffraction (XRD) measurements 
carried out at room temperature. The crystalline structure 
was refined using the Rietveld technique. The lattice con¬ 
stants, the type of impurities, and the respective volume 
fractions are reported in Tabled In particular, the variations 
of the unit-cell parameters as a function of the nominal flu¬ 
orine content are plotted in Fig. [l] We note two important 
features: (i) the a-axis length is slightly reduced with dop¬ 
ing. This effect is almost Ln-independent and has been as¬ 
cribed to the presence of small oxygen deficiencies in sam¬ 
ples grown at high pressure ! 7 ! 14 ’ 16 (ii) The c-axis length de¬ 
creases linearly with increasing x(F), too. In this case, the 
replacement of oxygen with the more electronegative flu¬ 
orine implies a reduction of the unit-cell volume as the F 
content increases. This result is similar to that observed 
also in other F-doped 1111 series (see, e.g., Ref. for the 
Ce case) and demonstrates the successful F-to-0 substitu¬ 
tion. 


B. Magnetometry measurements 

DC magnetization measurements were performed by 
means of a superconducting quantum interference de¬ 
vice (SQUID) magnetometer (Quantum Design) on all the 
tested samples. Both magnetization vs. temperature (from 
2 up to 300 K) and isothermal magnetization measure¬ 
ments vs. applied field, m(H), at selected temperatures 
were carried out. The linear (i.e. purely paramagnetic) be¬ 
havior of m[H ) at both low and high temperatures (not 
shown) demonstrates the absence of diluted ferromagnetic 
impurities. This was confirmed also by fiSR measure¬ 
ments at high temperature (see below), where no signifi¬ 
cant asymmetry damping could be observed within the 5 ps 
data-acquisition window used. 

In Fig. [2] the zero-field cooled dc magnetization mea¬ 
surements carried out at fi 0 H = 3 T are plotted as in¬ 
verse susceptibility vs. temperature. Numerical fits of 



x 


FIG. 1. (color online). Evolution of the a— and c—axis lengths 
of NdFeAsOj_ x F x at 300 K as a function of the nominal fluorine 
content. Continuous lines are only guides to the eye. 


j(T) data with a Curie-Weiss law restricted to the high- 
temperature regime allowed us to determine Xo, the resid¬ 
ual temperature-independent contribution to the suscepti¬ 
bility. Successively, a linear fit of the data with l/[j(T) — 
X 0 ] provided the Curie constant C, from which the Nd 3+ 
magnetic moment /i e[T (in the free-ion approximation) 
could be extracted. The obtained values lie in the range 
Peff = 3.5-3.8 /i B , i.e., very close to the value expected 
for the Nd 3+ free-ion magnetic moment (3.62 // B ). No 



T[K] 


FIG. 2. (color online). Temperature dependence of the inverse dc 
susceptibility in NdFeAsOj^F*. Data were recorded in zero-field 
cooled (ZFC) mode in an applied field p 0 H = 3 T. Lines represent 
fits using the l/[x(T) — Xoi model for T > 50K (for consistency 
among the normal and superconducting samples). 

evidence of Nd antiferromagnetic order was found down 
to 2 K, in agreement with previous neutron diffraction 
studies,which show that a possible AF ordering of Nd 
moments most likely occurs below 2 K. 

At low temperature, the samples with higher fluorine 
doping, x = 0.08 and 0.25, display a clear superconducting 
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transition, as shown by the low-field (1 mT) ZFC magneti¬ 
zation data in Fig. [3] and summarized in Table [TT] The super¬ 
conducting shielding fraction^ of about 80% is most likely 
a lower bound, since at low doping values the magnetic 
penetration depth X L becomes comparable to the grain size 
(micrometer scale), hence resulting in a reduced Meissner 
response. This result demonstrates the bulk character of 
the superconductivity in these samples. 



FIG. 3. (color online). Temperature dependence of the volume 
magnetic susceptibility in NdFeAsCh.-Tx f° r x = 0-08 and 0.25 in 
an applied field of 1 mT, following zero-field cooling. 


C. Muon-spin spectroscopy 

Muon-spin relaxation measurements were carried out at 
the GPS spectrometer (7tM3 beam line) at Paul Scherrer 
Institut (PSI), Villigen, Switzerland. 100% spin-polarized 
muons were implanted uniformly over the sample vol¬ 
ume and successively the relevant decay positrons were 
detected. Once muons thermalize at interstitial sites (al¬ 
most instantaneously and without loss of polarization), 
they act as sensitive probes, precessing in the local mag¬ 
netic field B M with a frequency = y M /(27r)B^, where 
y M /(27i) = 135.53 MHz/T is the muon gyromagnetic ratio. 
Both zero-field (ZF) and longitudinal-field (LF) pSR exper¬ 
iments were performed. Due to the absence of externally 
applied fields, ZF-pSR represents the best technique for in¬ 
vestigating the spontaneous magnetism and its evolution 
with fluorine doping. On the other hand, LF-pSR measure¬ 
ments can be successfully used to single out the dynamic or 
static character of the internal magnetic fields 3 

1. Zero-field. fiSR measurements 

Figure [4] shows the short-time ZF asymmetry (i.e. the 
muon-spin precession signal), A(t), for all the tested NdFe- 
As0 1 _ x F x samples at both low and high temperatures. At 
high temperatures, i.e. above the Neel temperature T N , the 
asymmetry signal is practically flat, with no oscillations and 
with a negligible decay of the initial polarization, mainly 
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FIG. 4. (color online). ZF-pSR asymmetry data for NdFeAsOj.*- 
F* samples, collected at high and low temperatures. Solid lines 
represent numerical fits by means of Eqs. (JTJ and ([2]). See text for 
details. 


due to the randomly-oriented nuclear magnetic moments. 
When the temperature is lowered below T N , the asymmetry 
becomes highly damped, reflecting the presence of a spon¬ 
taneous magnetic order whose character depends on x (F): 
at low fluorine content we observe well-defined asymme¬ 
try oscillations which, as x(F) increases, change quickly to 
highly damped ones. This behavior is indicative of a long- 
range (static) uniform magnetic order. For x = 0.08 the 
damping becomes so high that the oscillatory behavior dis¬ 
appears altogether, to be replaced by a fast decaying signal, 
typical of a wide distribution of fields attributed to short- 
range magnetic order. Finally, in the x = 0.25 case, we 
find no difference between the low- and high-temperature 
asymmetries, which implies the absence of a magnetically 
ordered phased To go beyond the above qualitative re¬ 
sults and achieve a better understanding, the various ZF- 
pSR time-domain data were fitted using the function: 


A ZF (t) 

^ 0 ) 


a^e 


2 f(rBnt)+a e 2 


+ a,e 


-X,t 


( 1 ) 


where a^ B are the two transverse asymmetry components, 
while a L represents a longitudinal component, each with 
respective decay coefficients a AB and X L . The transverse 
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TABLE I. Unit cell parameters and impurity fractions in NdFeAsOj.Tx samples as determined by Rietveld refinement of x-ray diffraction 
data taken at room temperature. 


x(F) 

a-axis (A) 

c-axis (A) 

Pure phase (%) 

FeAs (%) 

NdOF (%) 

NdAs (%) 

Other (%) 

0.03 

3.9669(2) 

8.5602(5) 

98(2) 

- 

- 

- 

<2 

0.05 

3.9676(2) 

8.5557(5) 

98(2) 

- 

- 

- 

<2 

0.08 

3.9617(4) 

8.5570(10) 

93(3) 

4.11(9) 

2.3(1) 

0.38(5) 

0.21(1) 

0.25 

3.9536(1) 

8.5394(5) 

84(1) 

- 

9.8(0.4) 

5.7(2) 

0.5(1) 


components refer to the two distinct muon implantation 
sites, normally expected for the 1111 family: the most pop¬ 
ulated site, named A, located next to the FeAs planes, and 
a secondary one, named B, close to the oxygen atoms in 
the LnO planes^ Since only ~ 15% of the total implanted 
muons are located in B sites, their contribution to the total 
asymmetry is too small to allow the extraction of a second 
precession frequency^ Therefore, we describe it only by a 
Gaussian decay (a®). The corresponding longitudinal com¬ 
ponents for muons stopped in both sites share similar decay 
rates and, hence, can be merged into a single term, a L . 

The analytical form of the term /(t) depends on the F 
content: for x = 0.03 and 0.05 the best fit was obtained 
with /(t) = the zeroth order Bessel function, 

generally attributed to incommensurate long-range ordered 
systems^ On the other hand, for x = 0.08, no coherent 
oscillations could be detected [/(t) = 1] and the time- 
dependent asymmetry was better fitted using the expres¬ 
sion: 


A ZF (t) 

^(0) 


r ate^+af 


B t] + 


a, e 




( 2 ) 


The transverse component a T (t) was fitted by a simple ex¬ 
ponential decay, whereas for the longitudinal one, a L (t), a 
stretched exponential was chosen in order to distinguish 
possible dynamic contributions (vide infra), with [3 the 
stretching coefficient. All the considered samples were 



FIG. 5. (color online). Temperature dependence of the magnetic 
volume fraction for samples with x(F) = 0.03, 0.05, and 0.08. 
Solid lines represent numerical fits by means of Eq. ©. 


polycrystalline, therefore, the randomly-oriented internal- 
field model was invoked to calculate the low-temperature 
magnetic volume fraction: when the whole sample shows 
static magnetic order, on average V 3 of the muons expe¬ 
rience a field parallel to their initial spin direction and do 
not precess (longitudinal component), while the remain¬ 
ing 2/3 precess around the local field (transverse compo¬ 
nent). This is exactly what we find at low temperature 
for the samples x = 0.03, 0.05, and 0.08 (see Fig. [4]). 
Therefore, all these samples can be considered fully mag¬ 
netically ordered at low temperature. To track the tem¬ 
perature dependence of the magnetic volume fraction, V m , 
the following procedure was used. From the longitudinal 
component, as extracted from Eq. 0, one can calculate 
V m {T ) = | (l — a L ) • 100%. This is the case of the sam¬ 
ple x = 0.08 and, for temperatures below ~ 10 K, also of 
the samples x = 0.03 and 0.05. For the latter two, V m at 
higher temperatures could be extracted using the expres¬ 
sion V m (T) = (l —a p ara} ' 100%, where a para represents 
the muon asymmetry related to paramagnetic (i.e., non¬ 
magnetic) species, as resulting from a transverse-field (TF) 
/iSR experiment, typically performed at 3 mT. The obtained 
V m (T) data for x = 0.03, 0.05, and 0.08 cases are shown 
in Fig. [5] To determine the average Neel temperatures and 
the corresponding transition widths, the V m (T) data were 
fitted using the following phenomenological function: 


v m m 


1 

2 


1 - erf 


t-tA 

^t n J 


(3) 


Here erf is the cumulative error function (typical of a nor¬ 
mal distribution). The fit results are summarized in Ta¬ 
ble [TT] In general, as the fluorine content x(F) increases 
there is both a gradual decrease of the average T N and a 
progressive broadening of the transition. 

Fits of low-temperature asymmetry data with Eq. (Q]) al¬ 
low us to determine also the internal field values at the 
lowest temperature, as probed by muons implanted in sites 
A and B. The dependence of these fields on fluorine doping 
x (F) is shown in Fig. [6j As the fluorine content increases, 
the average field value decreases and its Gaussian 

distribution broadens. The x = 0.08 case deserves partic¬ 
ular attention: A(l) does not show any coherent oscilla¬ 
tion and the internal field can be modeled by a broadened 
Lorentzian distribution of fields, whose values range from 
zero up to A AB /y M . 

In Fig.[7]we show the temperature behavior of the main 
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FIG. 6. (color online). Doping dependence of the internal field 
values as probed by implanted muons in sites A (filled symbols) 
and B (open symbols) at 5 K. For the low-doped samples, x = 0.03 
and 0.05, both the local-field value B A at site A and its width 
AB a = cr A /r u (plotted as error bar) could be determined, whereas 

[A H- 

for the site B the field distribution width AB B = <J b /Yij, is shown. 
For the x = 0.08 sample, only AB a - b = X AB /are displayed. The 
open square represents the internal field at site A for the undoped 
compound, as reported in Ref. [l^. The dotted and dashed curves 
are guides to the eye. 


parameters obtained by fitting the time-dependent asym¬ 
metry by Eq. ([2]). We notice that: (a) the longitudinal re¬ 
laxation rate has a peak at about 30 K; (b) the stretching 
coefficient shows a broad minimum centered again at 30 K; 
(c) the field width probed by the implanted muons at site 
A is almost constant up to the onset of the magnetic order, 
while the field probed at site B is very low (2-5 mT) and 
almost undetectable above 30 K (not shown). 


TABLE II. Magnetic properties of the NdFeAsCh-j-F* family, as de¬ 
termined from pSR and dc magnetometry measurements. The T N 
of the undoped compound is taken from Ref. fl2l. 


x(F) 

r N (K) 

at n (K) 

F m ag (%) 

T C (K) 

0 

135 




0.03 

84(1) 

11(1) 

>90 

- 

0.05 

40(1) 

10(1) 

>90 

- 

0.08 

40(3) 

25(4) 

>90 

20(2) 

0.25 

- 

- 

- 

51(1) 


2. Longitudinal-field [iSR measurements 

The sample x = 0.08 is the most interesting one, since, 
at low temperature, bulk superconductivity coexists with 
an ordered magnetic state. The absence of muon-spin pre¬ 
cession and a large damping of asymmetry at short times 
clearly indicate the presence of a wide distribution of fields. 
However, it is not obvious whether this distribution is due 
to static fields, or to strongly fluctuating (i.e., dynamic) 



T[K] 


FIG. 7. (color online), x = 0.08 T-dependence of the parame¬ 
ters extracted by fitting the time-domain /xSR asymmetry data by 
means of Eq. (0. The dashed curves are guides to the eye. 


magnetic moments. To disentangle these two possibilities, 
we carried out longitudinal-field (LF) /tSR measurements 
at 30 K. This temperature represents a good compromise, 
since it is above the superconducting critical temperature 
T c n 20 K (see Fig. [3]), yet sufficiently low to admit an 
ordered volume fraction of about 70% (see Fig. [5j. In an 
LF-decoupling experiment^ an external magnetic field By is 
applied along the initial muon-spin direction. If By is of the 
same order as or higher than the internal static fields, then 
it will significantly affect the muon polarization through a 
“spin-locking” effect. On the other hand, in case of strongly 
fluctuating internal fields, the effect of the external field is 
barely noticeable. 

In Fig. [8] the time-dependent asymmetry is shown for a 
selection of By values. The time-domain data were suc¬ 
cessfully fitted by means of Eq. (O, with the field depen¬ 
dence of the resulting best-fit parameters for the longitu¬ 
dinal component being shown in Fig. [9] Several interest¬ 
ing and unexpected features can be noticed: (i) From the 
T-dependent analysis we know that the static local field 
can be as high as ~ 70 mT [see Fig. Oc)]. Once this is 
quenched by an applied field of 75 mT, the longitudinal 
asymmetry is fully recovered only for fields above 0.2 T, (it) 
the longitudinal relaxation rate X L oc 1 /T 1 has a peak at 
about 30-50 mT; (iii) the exponent p reaches 1 only at 
high fields. 


















6 



FIG. 8. (color online). Time-domain LF-asymmetry data taken 
at T = 30 K in the x = 0.08 sample for different values of the 
longitudinal field. Solid lines represent numerical fits by means 
ofEq. ((2). 


III. DISCUSSION 

An understanding of the F-doping effects in Ndllll com¬ 
pounds, both regarding their magnetic and superconduct¬ 
ing properties, as well as the possible presence of an M-SC 
crossover region, is expected to provide a more complete 
picture of the iron-arsenide superconductors. The results of 
our systematic investigation of the structural and magnetic 
properties of representative compounds of the Nd family 
are summarized and discussed in detail below. 


A. Static magnetism and superconductivity 

In Ndllll the macroscopic magnetic behavior (dc sus¬ 
ceptibility) is dominated by the paramagnetic response of 
Nd 3+ ions, that display an AF order only at 1.96 K (in the 
undoped case) Contrary to the case of Cel 111,- where 
the AF order occurs at higher (i.e., more accessible) tem¬ 
peratures, the rather low values of the Nd family prevented 
us from verifying possible correlations between T N (Nd) and 
F-doping. Given these difficulties, we focused on the mag¬ 
netism of the Fe 2+ ions, crucial in determining most of the 
properties of the Lnllll compounds. 

Below T N (FeAs) (ca. 135 K), a commensurate, long- 
range AF order develops in the FeAs planes, giving rise to 
a coherent precession of the implanted muon ensemble^ 
However, already small values of F-doping are sufficient to 
induce a high damping in the oscillating muon polariza¬ 
tion. This is exactly the case of our x(F) = 0.03 and 0.05 
samples, where the oscillating term is well described by a 
zeroth-order Bessel function, normally the fingerprint of an 
incommensurate long-range magnetic orden^S^ when the 
F content is further increased, the size of the AF-ordered do¬ 
mains decreases and the magnetic damping at short times 
becomes so strong that no oscillations are detected in the 
time-dependent muon polarization. This situations is typ¬ 
ical of a short-range magnetic order— and this is exactly 



B lf [ mT ] 

FIG. 9. (color online). Longitudinal-field dependence of the pa¬ 
rameters extracted by fitting the time-domain asymmetry data 
taken at T = 30 K in the x = 0.08 sample (see Fig. [8j. The 
dashed curves are guides to the eye. 


what we found in the case of the sample with x(F) = 0.08. 

By considering the T N and T c values inferred by our de¬ 
magnetization and pSR data (see Table [TT]) one can draw 
a tentative phase diagram that describes the evolution of 
both M and SC orders as a function of F-doping. As shown 
in Fig. M the main feature of the phase diagram is the 
presence of a narrow region where electron doping induces 
both bulk superconductivity and FeAs magnetic order in 
the whole sample volume. The most reasonable way of 
reconciling the simultaneous presence of M and SC bulk 
phenomena is to assume their coexistence at a nanometer 
length scale, as already invoked in the case of the Smllllj^ 
Cel 111 — 12£ and Lai 111— families. Therefore, one can most 
likely conclude that this particular kind of M-SC coexis¬ 
tence is a common feature of all the Lnllll iron pnictides, 
independent of the lanthanide ion. 

On the other hand, contrary to the above cases, the 
x = 0.08 sample of the Ndllll family shows clear signs of 
non-negligible magnetic fluctuations near the FeAs ordering 
temperature. This peculiar feature is discussed in detail in 
the following section. 
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FIG. 10. (color online). Magnetic ordering T N and critical super¬ 
conducting temperatures T c vs. fluorine doping x in the NdFeAs- 
Oi-jcFjc family. The lines are guides to the eye. The T N value for 
the x = 0 case is taken from Ref. E 


B. Fluctuating magnetism in the x = 0.08 case 

The muon LF asymmetry data shown in Fig. [8] exhibit 
two clearly distinct relaxation regimes (see, e.g., data at 
B = OmT). At short times, the (fast) relaxation is domi¬ 
nated by static magnetic moments (on the /iSR time scale), 
while at long times, the (slow) asymmetry relaxation rate 
(oc 1 / Tj ) is mostly determined by dynamic (i.e., fluctuat¬ 
ing) magnetic moments. To account for both contributions, 
the LF-pSR data were fitted by means of Eq. (O, where the 
stretched exponential term can accommodate different sit¬ 
uations. Indeed, the stretching coefficient (3 can assume a 
value /3 = 2, indicating a Gaussian relaxation process (with 
the implanted muons probing quasi-static electronic mo¬ 
ments with short-range magnetic correlations), or /3 = 1, 
indicating a single spin-spin autocorrelation time for the en¬ 
tire electronic spin system. The case /? < 1, instead, is typ¬ 
ical of broad inhomogeneous distributions of dynamic (i.e., 
fluctuating) local fields (or relaxation rates), which can be 
due to randomness, disorder, or frustration, as e.g., in the 
case of spin glasses. 

Figures [Tfb) and (9fc), which report the evolution of the 
13 parameter with temperature and LF field, respectively, 
show similar features: [3 exhibits a broad minimum at 
about T N /2, where it assumes the value 0.6. Then, for tem¬ 
peratures below 20 K, where the magnetic moments in the 
FeAs planes become static, it grows back again to ca. 1. 
Also [3 = /3(B Lf ) shows a broad minimum centered in the 
0.05-0.1 T field range, where the longitudinal relaxation 
rate X L has a narrow peak reflecting the competition be¬ 
tween the static and dynamic components of the local field 
at the muon site. For applied fields above 75 mT, the in¬ 
dependence of the longitudinal relaxation rate is no more 
due to decoupling, but is determined only by the long-lived 
local field fluctuations. This fact is confirmed by the slow 
recovery of the longitudinal muon-asymmetry component. 
A similar behavior was observed also in Tb 2 Sn 2 0 7 ,— where 


the occurrence of such a peak was attributed to the inappli¬ 
cability of the Redfield relaxation model and, consequently, 
reported as evidence of a broad distribution of spin-spin 
correlation times. 

From these observations we argue that in the 20-40 K 
temperature interval the local magnetism shows a spin- 
glass character, most likely due to the proximity of the 
magnetic transition at T N ~ 40 K. Consequently, in this 
T range one can tentatively use the same approach valid 
for spin-glass systems. This will allow us to discuss the 
relation between the longitudinal-field data and the dy¬ 
namic spin fluctuations occurring in the magnetic phase. 
To this aim we recall that the relaxation at long times 
probes the spin-fluctuation dynamics and, in absence of 
cross-correlations (as is the case for spin glasses), the re¬ 
laxation rate is just the Fourier transform of the dynamic 
spin-spin autocorrelation function, q(t) = (S(0) • S(t)). 
When the spin-glass transition is approached from above, 
q(t) exhibits either a power-law (PL), q(t) oc t~ a , or a 
stretched-exponential (SE) dependence, q(t) oc e~^ Xt ^ .— 
Since in our case the magnetic transition is quite broad (see 
x = 0.08 data in Fig. O, we can use the same theoretical 
approach valid in spin glasses. The Fourier transform of 
the autocorrelation function, 5{q(t)}, leads then to a time- 
field scaling of the pSR asymmetry A(t,B LF ) = A(t/B[ F ), 
where y < 1 for power-law and y > 1 for stretched- 
exponential correlationsj2&22 This scaling stops being valid 
at high fields, where q(t) is directly affected by the applied 
LF field, typically for p B B LF > k B T , a limit which is never 
reached in our case. 

As a next step, one can try to establish the nature of the 
spin-spin correlation function and, from this, the nature of 
the dynamic component of the relaxation. By using the 
same method adopted in spin-glass systems, in Fig, fill we 
display the time-dependent asymmetry as a function of the 
scaled variable t/B[ F . To exclude from this analysis the ef¬ 
fects of the static component of the local field at the muon 
site, we removed the short-time data (i.e., t < 1 ps). With 
a value of y = 0.5 the t -dependent muon asymmetries ob¬ 
tained at different fields overlap over almost two orders 
of magnitude in t/B ] L b < 28 1 29 , except for very low field val¬ 
ues, where the static local field is expected to dominate. 
Since the scaling exponent y is lower than one, it sug¬ 
gests that at the pSR time scale the spin-spin correlation 
function q(t) is best approximated by a power law, rather 
than by a stretched exponential, hence suggesting a pic¬ 
ture of cooperative critical spin fluctuations rather than a 
distribution of local fluctuation rates. This finding seems in 
marked contrast with the field- and temperature-behavior 
of the stretched exponential parameter, [3 < 1, which hints 
at a broad inhomogeneous distribution of the local mag¬ 
netic fluctuation rates. The scenario of an inhomogeneous 
sample is, nevertheless, very unlikely: (i) the presence of 
a peak in A f (T) at about 30 K, that should be absent in a 
fully disordered system; (ii) at T = 30 K not less than 70% 
of the sample volume shows static ordered magnetism, and 
(iii) as suggested by the x-ray analysis and confirmed by 
a rather sharp superconducting transition (relative to this 
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FIG. 11. (color online). Muon asymmetry data taken at T = 30 K 
in the x = 0.08 sample in various applied longitudinal fields vs. 
the scaled variable t/B [ F , with y = 0.5 for x = 0.08. 


family), the x = 0.08 sample is expected to be fairly homo¬ 
geneous. Therefore, as in the case of CeFePO,— the sizable 
dynamic component evidenced in our sample most likely 
reflects a cooperative scenario. 

Interestingly, the critical slowing down of spin fluctuations 
occurs at sufficiently low temperatures, i.e., at T <20 K, 
where both the static magnetism and the superconductiv¬ 
ity occur over the whole sample volume. Unfortunately, we 
have no means of distinguishing whether the dynamic com¬ 
ponent can be related to the superconducting state or, most 
likely, it reflects a glassy magnetic phase. 

Finally, we remark that the disordered magnetic phase 
seems a constant feature of the M-SC crossover region of 
the iron-based superconductors that was observed, for in¬ 
stance, also in iron chalcogenides at intermediate Se-Te 
substitutions^ This type of short-range (disordered) mag¬ 
netism could be driven by orbital nematic correlations in¬ 
duced by the tendency towards ordering of the iron 3d 
orbitals^— To address this fundamental issue and unveil 


the true nature of the disordered magnetic phase, a sys¬ 
tematic study of the structural and of the in-plane transport 
properties vs. temperature on high quality NdFeAsO|_ x F.. 
single crystals should be seriously considered. 


IV CONCLUSION 

In this work we reported on the evolution with F-doping 
of the magnetic and superconducting properties of the 
Ndllll family. We have shown that its electronic phase di¬ 
agram presents similar features to other compounds of the 
Lnllll family (Ln = La, Sm, Ce). In particular, at low dop¬ 
ing (x < 0.06) a long-range antiferromagnetic order sets in 
for T < T n . At intermediate doping (x 0.08) a short- 
range magnetic order develops which, at low temperature, 
coexists with bulk superconductivity on a nanometer length 
scale. Finally, at high F-doping values (x 0.25), only the 
superconductivity survives, with no relevant traces of any 
magnetic order. 

Contrary to the other Lnllll cases, in the intermediate 
F-doping regime of Ndllll, we find clear evidence of fluc¬ 
tuating magnetism on the pSR time scale. This dynamic 
magnetic component persists down to low temperature, 
until a short-range static magnetism and superconductiv¬ 
ity develop over the whole sample volume. This peculiarity 
of Ndllll seems to be related to the glassy nature of its 
magnetically ordered phase at intermediate temperatures, 
i.e., just below the average Neel temperature. 
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